Electroweak Radiative Corrections to Muon Capture 



o 
o 

(N 



< 

o 



o 



X 



Andrzej Czarnecki 
Department of Physics, University of Alberta 
Edmonton, AB TOG 2J1, Canada 

William J. Marciano 
Brookhaven National Laboratory, Upton, NY 11973, USA 

Alberto Sirlin 

Department of Physics, New York University, 
4 Washington Place, New York, NY 10003, USA 



Electroweak radiative corrections to muon capture on nuclei are computed and found to be sizable. 
They enhance the capture rates for hydrogen and helium by 2.8% and 3.0% respectively. As a result, 
the value of the induced pseudoscalar coupling, g^ p , extracted from a recent hydrogen IS singlet 
capture experiment is increased by about 21% to g^ p — 7.3 ± 1.2 and brought into good agreement 
with the prediction of chiral perturbation theory, gp lcory = 8.2±0.2. Implications for helium capture 
rate predictions are also discussed. 
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The study of muon capture by nuclei, fi~N — > v^N , has played an important role in the development of weak 
' interaction physics [HQ. Used primarily in the past to explore nuclear structure and its effects on weak interactions, 
muon capture can now be employed to test Quantum Chromodynamics (QCD) and its basic chiral symmetries In 
addition, it can provide a possible window or constraint on new high mass scale physics Q, beyond Standard Model 
expectations, such as additional gauge bosons, charged Higgs scalars, leptoquarks etc. Of course, to be competitive 
i-C ■ with other precision low energy experimental tests of the Standard Model, both theory and experiment for muon 

capture must be known to a fraction of a percent. 
H ■ Here, we would like to advance the theory of muon capture to that high level of precision by including Standard 
J> ' Model electroweak radiative corrections and estimating their degree of reliability. From our previous work [a, El 0, H, HI 
on neutron (and nuclear) /3-decay, one can anticipate that such quantum loop effects are relatively large, ~ 2 — 3%, 
and therefore important for any precision confrontation between muon capture theory and experiment. As we shall 
show, that indeed is the case. 

We begin by recalling the basics of muon capture. Negative muons, /i - , are stopped in matter. They bind 
electromagnetically with nuclei and quickly cascade down to the lowest energy atomic orbitals. There, primarily 
from 15 states, the muon's final fate is to undergo either ordinary muon decay, fi~~ — > e~i? e u^, or weak capture 
fi~ N — » v^N' on the nucleus. 

Ordinary decay in orbit occurs essentially at the same rate as in vacuum (modulo bound state time dilation and 
other small effects 

OH G3). The already well known "free" muon lifetime has been recently remeasured [12j, thereby 
leading to the improved world average 



r p = 2.197019(21) x 1CT 6 sec (1) 

Further improvement by an additional factor of ten is expected. 

The competing weak capture reaction, pT N — > v^N 1 ', proceeds via W boson exchange with the nucleus. Due 
to an overlap flux factor from its atomic wavefunction at the origin squared and a factor of Z (nuclear charge) 
corresponding to the number of nuclear protons that can induce capture, the overall capture rate scales very roughly 
as Z . In hydrogen [Z = 1), the capture rate is predicted to be very small. From the 15 singlet (spin 0) /j,~p 
state, it is only about 0.16% of the ordinary decay rate and for the triplet (spin 1) bound state configuration, it is 
a tiny 0.0025%. Those small rates make experimental hydrogen capture studies difficult, which is unfortunate, since 
hydrogen theory is very clean. Decay and capture rates become comparable for Z ~ 10, while at much higher Z , 
capture dominates. 

An interesting technique used to obtain muon capture rates involves comparing free and bound fi~ lifetimes, 

r (M"^V - v»N>) = -jJLj - -1^, (2) 



(after making small bound state lifetime corrections). Using an ingenious application of that lifetime technique, the 
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MuCap collaboration [13| at PSI recently reported a precise measurement of the IS singlet capture rate in hydrogen, 

r(jjTp -> i/^) 8 ;*' = 725.0 ± 13.7(stat) ± 10.7(sys)/sec. (3) 

That already impressive ±2.4% level of accuracy is expected to further improve to better than ±1% as additional 
data is analyzed. 

In the case of helium, the capture rate for fj," 3 He — ► 3 H has been even better measured by directly detecting the 
charged final state 3 H. For the statistical combination of singlet (spin 0) and triplet (spin 1) 15 pT 3 He bound states, 

T(p- 3 Hc -» i/„ 3 H) stat = ir(/i" 3 He -» ^ 3 H) singlct + -1^" 3 He - ^ 3 H) triplet , (4) 

a long standing result [13] 

r(/i" 3 He -» ^ 8 H)£t = 1496(4)/ sec, (5) 

represents a remarkable ±0.3% determination. 

The Standard Model theoretical prediction for the basic p~p — » i/^n capture rate depends on four relativistic form 
factors that result from nucleon matrix elements of the V — A weak quark charged current, 



(n\dj a (l -~f 5 )u\p) = u n {p 2 ) 
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(6) 



Two other form factors, scalar and pseudotcnsor, arc in general possible, but are negligibly small in the Standard 
Model (arising from isospin violation). They should, however, be included in general searches for "New Physics" 
effects [3,[l5j]. In terms of the above form factors, the capture rate is given by (modulo radiative corrections, discussed 
later) 



r(/i p -> Vpri) 



singlet 



T(p P^^n)| triplet 
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(7) 



M denotes the mass of the p p atom. We neglect the binding effect and use M = m p + nip. The p p hydrogenic 
wave function at the origin is 



|^(0)| 2 = (1 - Aapr p ) ~ (1 - 0.005) , 



7T 

m p nip 
m p + m p 



(reduced mass) 



(8) 



where we have accounted for the proton charge distribution with the radius r p = °^j 2 fm (see (l6l . [l7T | for a more 
detailed discussion). 

Three of the four form factors in Eq. ^ are very well determined at q 2 = from CVC and neutron (3 decay [1, [l8| , 



Fi(0) = 1, 
F M (0) = 3.706, 
3 a(0) = 1.2695(29). 



(9) 
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Extrapolating to q 2 = — 0.88m 2 , as appropriate for /i capture on hydrogen, one finds 

F 1 (q 2 ) = 0.976(1), 
F M (q 2 ) = 3.583(3), 

g A (q 2 ) = 1.247(4), (10) 

where the errors include estimated q 2 evolution uncertainties. 

In the case of the induced pseudoscalar coupling, gp(qa), PCAC (partially conserved axial current) and chiral 
perturbation theory predict 0, [H HH, H2| 

2 2m M <7 7rp „ ((7q )Fn ^ fn\ 2 f~ l ~ l \ 

9p{%) = 5—2 ■z9A{0)m IJ .mNr A , (11) 

which for 57rpn = 13.05(20), = 92.4(4) MeV, and r\ = 0.43(3) fm 2 implies 

g P (q 2 ) =8.2 ±0.2. (12) 

That prediction is expected to be very reliable, depending only on the chiral properties of QCD and principles of 
PCAC. Nevertheless, it would be very useful to have a first-principles lattice QCD calculation of gp(q 2 )) (as well as 
9a{Qo))- Of course, it is also very important to verify the prediction in Eq. (|12p experimentally. 

Employing the above form factors at g 2 , and allowing for the variation gp{q$) = 8.2 + Sgp, one obtains from Eq. ([7]) 
the singlet 15 capture rate on hydrogen, 

T(p-p -> ^n)i5 Slct = 692.3(3.4) (1 + RC(H)) (1 - 0.0108<5 5P ) 2 / sec . (13) 

Gf, = 1.166371(6) x 10~ 5 GeV~ 2 (the Fermi constant obtained from the free muon lifetime [3), V ud = 0.9738 and a 
0.5% reduction from the finite proton size have been incorporated into Eq. (fl~3|) . The 1 + RC(H) factor represents the 
effect of electroweak radiative corrections, which up until this work have not been seriously considered in discussions 
of muon capture [23|, [24[. If we set RC(H) = and compare Eq. (p~3|) with Eq. ([3]), we find gp(q 2 ) = 6.0 ± 1.2 which is 
about 2a below the prediction in Eq. (|12p : however, that result is not very meaningful since we expect the radiative 
corrections to be sizable. 

In the case of helium, the tree level theoretical prediction for muon capture is not as pristine. When compared 
with the same input parameters, two distinct approaches give somewhat different results. The first is based on an 
elementary particle prescription which treats 3 He and 3 H as initial and final particle states [!, [25|. It then employs 
form factors analogous to those in Eq. © (but defined with an additional minus sign for all but Fi) at q 2 = —0.954m 2 
appropriate for \i~ capture on 3 He — > 3 H. Using CVC for the vector form factors and PCAC to relate axial-vector and 
pseudoscalar form factors, the analysis leads to what has been viewed as a rather reliable 3 He capture rate prediction. 
It depends primarily on the input 

g A (q 2 = -0.954m 2 )3 Hc ^3 H - 1-052 ± 0.005, (14) 

obtained by evolving <7a(0)3h o ^ 3 h = 1-212, obtained from tritium {3 decay [26|,[23], to q 2 = —0.954m 2 . 

The second method for calculating the capture rate for /i~ 3 He — * 3 H uses an impulse approximation to combine 
the basic fi~p — + v^n captures within 3 He p, 28]. It has been argued that when supplemented by meson exchange 
current corrections [2^|, this method agrees with the above (elementary particle) approach. However, a close scrutiny 
of the most detailed impulse approximation study [30| reveals some difference in their predictions. 

Normalizing to V u a — 0.9738, the elementary particle model approach predicts [28| 

T ( M - 3 Hc -> ^ 3 H) s E ™ = 1492(21) • (1 + RC(He))/ sec, 

while the impulse approximation study by Marcucci et al. [3(| updated to a central value of g A 
g P = 8.2 + 5g P gives 

r (aT 3 He -> i/„ 3 H)^ at = 1462(8)(7) 3A • (1 + RC(He))(l - 0.013<5. 9P )/ sec, (16) 

Again, we have allowed for inclusion of electroweak radiative corrections, RC(He), appropriate for capture. For clarity, 
we note that the values reported by Marcucci et al. [3(| are larger than in Eq. (fTB| because these authors identified Gy 
with a parameter Gy = 1.024|l/ u d| 2 G 2 , extracted from superallowed beta decays, in which inner radiative corrections 



(15) 

= 1.2695(29) and 
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of 2.4% were already included. In Eq. (|16p . we have factored out this 2.4% effect and included it in the overall RC(He) 
to be discussed below. 

The prediction in Eq. (fT"5|) is in very good agreement with Eq. ([5]), T(fi~ 3 He — > f M 3 H)g£f t = 1496(4)/ sec, if we 
naively set RC(He) = 0. That agreement has been viewed as a success of theory and used to constrain [4J "New 
Physics" appendages to the Standard Model. On the other hand, Eq. (fT6|) only agrees with experiment if one includes 
the +2.4% radiative correction contained in their G\ value. 

Now, we consider the electroweak radiative corrections (RC). They naturally divide into two contributions. The 
first set is essentially common to all semileptonic weak charged current amplitudes normalized in terms of G^, the 
Fermi constant obtained from the free muon lifetime. The second type of correction stems from QED corrections 
to the muonic atom wavefunction. As pointed out by Goldman (3lj , those latter effects are dominated by vacuum 
polarization corrections to the Coulombic bound state interaction. 

Making the above division, 

RC(JV) = RC(iV)i + RC(JV) 2 , (17) 

we find from the detailed studies of neutron decay [1, Q (neglecting terms of relative order am M /mjv) that the 0(a) 
electroweak radiative corrections to the muon capture rate on hydrogen are given by 



RC(H)i = f 

Z7T 

where m z = 91.1875 GeV, m p = 0.938 GeV, 



4 In — - 0.595 + 2C + 5 (tou,/3u = 0) 



(18) 



C = 0.829, (19) 



and the quantity g{m tl ,(3 ll — 0) can be obtained from Eq. (20b) in Ref. |32| by replacing m e — * m^, ignoring 
bremsstrahlung and taking the non-relativistic (zero muon velocity) [3^ = limit. In that way one finds 

<?(m M ,/3 M = 0) =3hi^ = -0.199. (20) 
m M 4 

In total, Eq. (fl~8|) gives 0.0223. Summing up higher order leading logs along the lines of ref. @, Q enhances that 
correction somewhat to 

RC(H)i = 0.024(4), (21) 

where we have included a fairly generous estimate of the uncertainty. It corresponds to roughly a ±100% variation in 
C and conservatively allows for 0{am )Jj /mp) corrections that we have not computed. 

We note that the first two bracketed terms in Eq. fTg)) (which include QCD perturbative effects) are of short-distance 
origin and therefore apply to all muon capture rates. Similarly the g function is essentially unchanged as long as the 
muon is non-relativistic and 0{am^/mp) contributions are ignored. On the other hand, the quantity C in Eq. (|19p 
is specific to hydrogen and will be modified by nucleon interactions in multi-nucleon systems. Rather than try to 
account for that modification, we assume that our rather conservative error covers those variations and continues to 
hold, 

RC(He)i = 0.024(4), (22) 

If needed, the correction in Eq. (|22[) can be used as a good approximation for any muon capture rate. We note that 
our +2.4% correction happens to coincide numerically with the inner radiative corrections included in the G v value 
employed in [30l ]. 

At this point we note that the factorization of the radiative corrections comes about because in the formulation 
of Ref. 8] , which we follow, the axial couplings in neutron decay have by definition the same electroweak radiative 
corrections as the vector ones. Small differences that can result from q 2 + 1 are included in the theoretical uncertainty 
or evolution uncertainty of the form factors. 

The vacuum polarization correction to the muon bound state wavefunction [3l| must be individually evaluated for 
different nuclei. A detailed calculation gives 

RC(H) 2 = 1.73- ~ 0.004, (23) 
which is somewhat smaller than found by Goldman [3l| . In the case of helium, we obtain 

RC(He) 2 = 2.92- ~ 0.0068. (24) 

7T 
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Overall, we find 



1 + RC(H) = 1.028(4), 

l + RC(He) = 1.030(4), (25) 



which modify the capture rate predictions in Eqs. (fl~3|) and (|15H16I) to 



r(M"P -> ^")is Slet = 711.5(3.5) SA (3) flC (l-0.01085 5 p)7se C) 



. EPM 



r(/x- 3 He^^ 3 H);;;" = 1537(22)/ sec, 
r(/x- 3 He^^ 3 H)^ at = 1506(8)(7) gA (6) RC (l-0.0135 3P )/sec. (26) 

For hydrogen, comparison with the experimental results in Eq. leads to 5g p = —0.9 ± 1.2, 

g c p p = 7.3 ±1.2 hydrogen. (27) 

The electroweak radiative corrections have increased the value of g e p P by about +21%. They bring theory and 
experiment into agreement. That situation is to be contrasted with the world average g C p P — 10.5 ± 1.8 obtained [|| 
from muon capture on hydrogen before the new MuCap result fl3T | and our evaluation of the radiative corrections. 
On its own, that previous world average would have been shifted to g C p P — 11.7 ± 1.8 by the radiative corrections, 
about a 2er deviation from the chiral perturbation theory prediction. However, including the MuCap fl3j result, one 
finds the new world average from muon capture on hydrogen, gp* p = 8.7 ± 1.0, in good agreement with the theoretical 
prediction of Eq (fT2|) . 

For helium, radiative corrections spoil somewhat the good agreement between experiment and the EPM prediction. 
The new disagreement suggests a smaller value of <m(~ 0.954m^)H C ^H is likely or a significantly larger gp by about 
25% in magnitude beyond PCAC predictions (a situation similar to hydrogen if we had used the pre MuCap capture 
rates). On the other hand the IA approach [3(| fares much better, leading to g C p P — 8.7 ± 0.6 which is also in good 
agreement with chiral perturbation theory. 

In summary, when our calculation of the electroweak radiative corrections to muon capture on hydrogen is combined 
with a new singlet /i~p capture rate measurement, it leads to g C p P — 7.3 ±1.2 which is in very good accord with the 
prediction of chiral perturbation theory, <7p Cory = 8.2 ± 0.2. That agreement would seem to close a confusing chapter 
in nuclear physics which has seen decades of disagreement regarding the value of g°p P ■ It will be very interesting to 
watch continuing improvements in the MuCap results. 
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